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Abstract —Attitude estimation is often a prerequisite for 
control of the attitude or orientation of mechanical systems. 
Current attitude estimation algorithms use coordinate rep¬ 
resentations for the group of rigid hody orientations. All 
coordinate representations of the group of orientations have 
associated problems. While minimal coordinate representations 
exhibit kinematic singularities for large rotations, non-minimal 
coordinates like quaternions require satisfaction of extra con¬ 
straints. A deterministic attitude estimation problem for a rigid 
body with bounded measurement errors is considered here. 
An attitude estimation algorithm that globally minimizes the 
attitide estimation error, is obtained. Assuming that the initial 
attitude, the initial angular velocity and measurement noise 
lie within given ellipsoidal bounds, an uncertainty ellipsoid 
that bounds the attitude and the angular velocity of the 
rigid body is obtained. The center of the uncertainty ellipsoid 
provides point estimates, and the size of the uncertainty ellipsoid 
measures the accuracy of the estimates. The point estimates, and 
the uncertainty ellipsoids are propagated using a Lie group 
variational integrator, and its linearization, respectively. The 
attitude estimation is optimal in the sense that the attitude 
estimation error and the size of the uncertainty ellipsoid is 
minimized. 

I. Introduction 

Attitude estimation is often a prerequisite for controlling 
aerospace and underwater vehicles, mobile robots, and other 
mechanical systems moving in space. Hence, attitude esti¬ 
mation may be used in spacecraft and aircraft, unmanned 
vehicles and robots, including walking robots. In this paper, 
we look at the attitude estimation problem for the uncon¬ 
trolled dynamics of a rigid body in an attitude-dependent 
force potential. The estimation scheme we present has the 
following important features; (1) the attitude is globally 
represented without using any coordinate system, (2) the 
filter obtained is not a Kalman or extended Kalman filter, 
and (3) the attitude and angular velocity measurement errors 
are assumed to be bounded, with ellipsoidal uncertainty 
bounds. The static attitude estimation using a global attitude 
representation is based on [1]. Such a global representation 
has been recently used for partial attitude estimation with a 
linear dynamics model in [2]. 

The attitude determination problem for a rigid body from 
vector measurements was first posed in [3]. A sample of 
the literature in spacecraft attitude estimation can be found 
in [4], [5], [6], [7], [8]. Applications of attitude estimation 
to unmanned vehicles and robots can be found in [2], [9], 
[10], [11]. Most existing attitude estimation schemes use 
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coordinate representations of the attitude. As is well known, 
minimal coordinate representations of the rotation group, like 
Euler angles, Rodrigues parameters, and modified Rodrigues 
parameters (see [12]), usually lead to geometric or kinematic 
singularities. Non-minimal coordinate representations, like 
the quaternions used in the quaternion estimation (QUEST) 
algorithm and its several variants ([4], [8], [13]), have their 
own associated problems. Besides the extra constraint of 
unit norm that one needs to impose on the quaternion, the 
quaternion representation for a given rotation depends on the 
sense of rotation used to define the principal angle, and hence 
can be defined in one of two ways. 

A brief outline of this paper is given here. In Section II, the 
attitude determination problem for vector measurements with 
measurement noise is introduced, and a global attitude deter¬ 
mination algorithm which minimizes the attitude estimation 
error is presented. In Section III, the attitude dynamics and 
dynamic estimation problem is formulated, and an algorithm 
to numerically integrate the dynamics is presented. Section 
IV presents the attitude estimation scheme with attitude and 
angular velocity measurements. Section V presents some 
simulation results followed by conclusions in Section VI. 

II. Attitude Determination 
A. Attitude determination from vector observations 

Attitude of a rigid body is defined by the orientation of 
a body fixed frame with respect to a reference frame, and 
the attitude is represented by a rotation matrix that is a 
3x3 orthogonal matrix with determinant of 1. Rotation 
matrices have a group structure denoted by SO (3). The group 
operation of SO(3) is matrix multiplication, and its action 
on takes a vector represented in body fixed frame into 
the reference frame by matrix multiplication. 

We denote the known direction vector of the ith point 
in the reference frame as e* G S^, and the corresponding 
vector represented in the body fixed frame as 6* G S^. 
These direction vectors are normalized so that they have 
unit lengths. The e® and V are related by a rotation matrix 
C G SO(3) that defines the rigid body attitude; 

e* = Cb\ 

for all i G {1,2,-,to}, where m is the number of 
measurements. We assume that e® is known accurately and 
6® is measured by sensors in the body fixed frame. Let the 


measured direction vectors (with sensor errors) be denoted 
6 * G and let an estimate of the rotation matrix be denoted 
C G SO(3). The estimation error is given by 

e* - db\ 

The attitude determination problem consists of finding an 
estimate C G SO (3), and is given by the following weighted 
least squares problem: 

^ m 

minj"=-y^ Wi(e* — C 6 *)'^(e* — (76*), 

= itr {E-CBfw{E-CB) , (1) 

subject to (7 G SO(3), 

where E = [e^, e^, • • • , e"*] G 

B = [61,62,-•• , 6 ™] G and W = 

diag [tt;!, , ty™] g has weight factors 

for each measured vector. 

This problem is known as Wahba’s problem [3]. The 
solution in terms of quaternions, known as the QUEST 
algorithm, is presented in [7]. A solution without using 
generalized attitude coordinates is given in [1]. A necessary 
condition for optimality of ( 1 ) is given by 

^ ( 2 ) 

where L = EWB'^ G 

The following result, which is proved in [1], gives an 
unique estimate C G SO(3) of the attitude matrix that solves 
the attitude determination problem ( 1 ). 

Theorem 1: The unique minimizing solution to the atti¬ 
tude determination problem ( 1 ) is given by 

C = SL, S = (3) 

where 

L = QR, Q G SO(3), (4) 

and R is upper triangular and invertible; this is the QR de¬ 
composition of L. The symmetric positive definite (principal) 
square root is used in (3). 

The proof is based on the fact that is a Morse function, i.e., 
its critical points are non-degenerate. From the Morse lemma 
[14], we conclude that these non-degenerate critical points 
are isolated, and hence the estimate given by (3) uniquely 
minimizes the attitude estimation error. 

B. Estimation with bounded state uncertainties 

A stochastic state estimator requires probabilistic models 
for the state uncertainty and the noise, which are often not 
available. Assumptions are usually made on the statistics 
of disturbance and noise processes, in order to make the 
estimation problem mathematically tractable. In many prac¬ 
tical situations such idealized assumptions are not appro¬ 
priate, and may cause poor estimation performance [15]. 
An alternative deterministic approach is to specify bounds 
on the uncertainty and the measurement noise without any 


assumptions on their distribution. Noise bounds are available 
in many cases, and such a deterministic estimation scheme 
is robust to the noise distribution. An efficient but flexible 
way to describe the bounds is using ellipsoidal sets, referred 
to as uncertainty ellipsoids. 

This deterministic estimation procedure for a 2 dimen¬ 
sional system is illustrated in Fig. 1, where the left figure 
shows time evolution of an uncertainty ellipsoid, and the 
right figure shows a cross section at a fixed time when the 
state is measured. At the fcth time step, the state is bounded 
by an uncertainty ellipsoid centered at Xk- This initial ellip¬ 
soid evolves over time. Depending on the dynamics of the 
system, the size and the shape of the tube are changed. At 
the k -f 1 th time step, the predicted uncertainty ellipsoid is 
centered at xl^^- The state is then measured by sensors, 
and another ellipsoidal bound on the state is obtained by 
the measurements. The measured uncertainty ellipsoid is 
centered at The state lies in the intersection of the 

two ellipsoids. In the estimation procedure, we find a new 
ellipsoid that contains the intersection, which is shown in 
the right figure. The center of the new ellipsoid, Xk+i is 
considered as a point estimate at time step k + 1, and the 
magnitude of the new uncertainty ellipsoid measures the 
accuracy of the estimation. This deterministic estimation is 
optimal in the sense that the size of the new ellipsoid is 
minimized. 

A deterministic estimation process based on set theoretic 
results was developed in [16]. Optimal deterministic esti¬ 
mation is considered in [17] and [18], where an analytic 
solution for the minimum ellipsoid that contains a union or 
an intersection of ellipsoids is obtained. 

III. Attitude Dynamics and Dynamic Attitude 
Estimation 

A. Equations of motion 

We now consider dynamic state estimation of the attitude 
dynamics of a rigid body in a potential U(C) : SO (3) R 
determined by the attitude, (7 G SO(3). A spacecraft on a 
circular orbit including gravity gradient effects [19], or a 3D 
pendulum [20] can be so modeled. The continuous equations 
of motion are given by 

JcIj -f w X Ju! = M, (5) 

(7 = CS{uj), (6) 

where J G R^^^ is the moment of inertia matrix of the rigid 
body, a; G R^ is the angular velocity of the body expressed 
in the body fixed frame, and S')-) : R^ so(3) is a skew 
mapping defined such that S{x)y = x x y for all x,y G R^. 
M G R^ is the moment due to the potential. The moment is 
determined by S{M) = ^ R — or more explicitly, 

M = n X Vr^ + r2 X Vr2 + ^3 X Vr ^, (7) 

where r*, G R^^^ are the ith row vectors of (7 and 
respectively. The derivation of the above equations can be 
found in [ 20 ]. 
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(a) Propagation of uncertainty ellipsoid 


(b) Filtering procedure 


Fig. 1. Uncertainty ellipsoids 


General numerical integration methods, including the pop¬ 
ular Runge-Kutta schemes, typically preserve neither first 
integrals nor the characteristics of the configuration space, 
SO(3). In particular, the orthogonal structure of the rotation 
matrices is not preserved numerically. To resolve these 
problems, a Lie group variational integrator for the attitude 
dynamics of a rigid body is proposed in [20]. This Lie 
group variational integrator is described by the discrete time 
equations. 

hS{Juk + ^Mk) = FkJd - JdFj, ( 8 ) 

Ck+i = CkFk, (9) 

J^^k+i = f'^ Jloi- + —F'^M^ + —( 10 ) 

where Jd S is a nonstandard moment of inertia matrix 
defined by Jd, — 5 tr[J]/ 3 x 3 — J, and F^ G SO(3) is 
the relative attitude over an integration step. The constant 
G K. is the integration step size, and the subscript k 
denotes the fcth integration step. This integrator yields a 
map {Ck,ujk) ^ (Cfc+ijCUfe+i) by solving (8) to obtain 
Fk G SO(3) and substituting it into (9) and (10) to obtain 
Ck+i and ujk+i- 

Since this integrator does not use a local parameteriza¬ 
tion, the attitude is defined globally without singularities. 
It preserves the orthogonal structure of SO(3) because the 
rotation matrix is updated by a multiplication of two ro¬ 
tation matrices in (9). This integrator is obtained from a 
discrete variational principle, and it exhibits the characteristic 
symplectic and momentum preservation properties, and good 
energy behavior characteristic of variational integrators. We 
use (8), (9), and (10) in the following development of the 
attitude estimator. 


uncertainty ellipsoid is a submanifold of TSO(3). An un¬ 
certainty ellipsoid centered at (C,a)) is induced from an 
uncertainty ellipsoid in R®, using the Lie algebra so(3); 


f(C,cD,P) 

= jc G SO(3), w G 

= jc'GSO(3), wGR^ 


U, 5lJ^ 


'C' 

5lo 


P- 


'C 

5uj 


< 1 


,P) 


( 11 ) 


where S{() = logm G so(3), Suj = uj — il) G R^, 

and P G is a symmetric positive definite matrix. 

Equivalently, an element {C,uj) G S{C,U!, P) can be written 
as 


C = Ce^^<\ 
OJ = (h + Soj, 


for X = 



G R® satisfying x^P ^x < 1. 


C. Measurement error model 

We give the measurement error models for the direction 
vector and for the angular velocity. The direction vector 
P G is in the body fixed frame, and let P G denote the 
cotTesponding measured directions. Since we only measure 
directions, it is inappropriate to express the measurement 
error by a vector difference. Instead, we model it by rotation 
of the measured direction; 

P = 

-P + S{F)b\ ( 12 ) 


B. Uncertainty Ellipsoid 

The configuration space of the attitude dynamics is SO (3), 
so the state evolves in TSO(3). Thus the corresponding 


where F G R^ is the sensor error, which represents the 
Euler axis of rotation vector from P to P, and ||i/*|| is the 
corresponding rotation angle in radians. The second equality 












assumes small measurement errors. The angular velocity 
measurement errors are modeled as 

Wfc = (ifc + Vk, (13) 

where Cjk S is the measured angular velocity, and Vk S 
is an additive error. 

We assume that the initial conditions and the sensor noise 
are bounded by prescribed uncertainty ellipsoids. 


(Co, Wo) G £{Co,ojo, Po), 

(14) 

F,€£^.i0,Sl), 

(15) 

fk G ^R3(0, Tk), 

(16) 


where Pq S Sl,Tk G are symmetric positive 

definite matrices that define the shape and the size of the 
uncertainty ellipsoids. 

IV. Attitude Estimation with Angular Velocity 
Sensor 

In this section, we develop a deterministic estimator for 
the attitude and the angular velocity of a rigid body as¬ 
suming that both attitude measurement and angular velocity 
measurements are available. The estimator consists of three 
stages; flow propagation, measurement, and filtered update. 
The propagation is to predict the uncertainty ellipsoid in the 
future. The measurement is to find an uncertainty ellipsoid 
in the state space using the measurements and the measure¬ 
ment error model. The filtered update finds a new estimate 
using the predicted uncertainty ellipsoid and the measured 
uncertainty ellipsoid. 

The subscript k denotes the fc-th discrete index. This may 
not coincide with measurement instants as we may resolve 
the evolution of the trajectory more frequently than the 
frequency of the measurements. This enables us to deal with 
measurements that are rather infrequent, with nontrivial at¬ 
titude evolution between the measurements. The superscript 
/ denotes the variables related to the flow update, and the 
superscript m denotes the variables related to the measure¬ 
ment update. ~ denotes a variable measured by sensors, and 
• denotes an estimated variable. 

A. Flow propagation 

Suppose that the attitude and the angular momentum at 
the kih step lie in a given uncertainty ellipsoid: 

{Ck ; G S{Ck j ^k^ ^ 

and suppose that new measurements are taken at the k -f /th 
time step. 

The flow update obtains the the uncertainty ellipsoid at the 
k -I- /th step using the given uncertainty ellipsoid at the fcth 
step. We assume that the given uncertainty ellipsoid at the 
fcth step is sufficiently small that the states in the uncertainty 
ellipsoid can be approximated by linearized equations of 
motion. This guarantees that the boundary of the state 
uncertainties at the k -f (th step remains an ellipsoid. 


Center: For the given center, {Ck,Cjk), the center of 
the uncertainty ellipsoid is obtained from the 

discrete equations of motion, (8), (9), and (10): 

hS{JuJk + ^Mfc) = FkJd - JdPj, (17) 

= CkPk, (18) 

J'^k+l ~ Pk Wfc -f '^Pk Mk + '^Mk+1- (19) 

This integrator yields a map (Ck,uJk) and 

this process can be repeated to find the center at the k -f /th 
step, 

Uncertainty matrix: The uncertainty matrix is obtained by 
linearizing the above discrete equations of motion. At the 
{k + l)th step, the state is given by perturbations from the 
center as 

Ck+i = 

uJk+i = 

for some C^+i: ^^l+i ^ Assume that the uncertainty el¬ 
lipsoid at the fcth step is sufficiently small. Then, 
are given by the following linear equations in [19]: 


-1 

r 


'4 

4 ' 

Ck 

[^4+i_ 


Ul 

41 

_5ujk_ 


where G can be suitably defined. 

Equivalently, we rewrite the above equation as 

4+1 = 4^k, 

where Xk = [Ck’^‘^k4 ^ 4 ^ Since 

{Ck,u!k) & £(Ck,Lbk,Pk), Xk G£R6(0,Pfc)by the definition 
of the uncertainty ellipsoid given in (11), we can show that 

4xk G ^0, ^ . 

Thus, the uncertainty matrix at the k + 1th step is given by 

Pi+i=AiPk[Aiy. (20) 

The above equation can be applied repeatedly to find the 
uncertainty matrix at the k + (th step. In summary, the 
uncertainty ellipsoid at the (k -I- 0th step is computed using 
(17), (18), (19), and (20) as: 

{Ck+i,ujk+i) G 

k+l ’ 4+1 ), (21) 

B. Measurement update 

The measured attitude and angular velocity have uncertain¬ 
ties due to sensor errors. However, we can find a uncertainty 
bound on the states because we assume that the sensor 
errors are bounded by known uncertainty ellipsoids. The 
measurement update obtains an uncertainty ellipsoid in the 
state space using the measurements and the sensor error 
models. 

Center: The center of the uncertainty ellipsoid, 
is obtained from the measurements. The 



attitude is determined by measuring the directions 
to the known points in the inertial frame. Let 
the measured directions to the known points be 
Bk+i = G Then, the attitude 

satisfies the following necessary and sufficient 
condition given in ( 2 ). 

(cr+y Lk+i - Ll+iCr+i = 0, ( 22 ) 

where Lk+i = Ek+iWk+iBj^i G The solution of 

(22) is obtained by a QR factorization of L^+i as given in 
Theorem 1. 

cr+z = Lk+u (23) 

where g G SO(3) and R G is upper triangular such 
that Lk+i = QR. The angular velocity is measured directly 
by sensors; 

QT+i=^k+i. (24) 

Uncertainty matrix: We can represent the actual state at 
the k + (th step as follows: 

Cfc+z (25) 

Wfc+i = (26) 

for Ck+h S The uncertainty matrix is obtained by 

finding an ellipsoidal bound for 

For the attitude, we transform the uncertainties in the 
directional sensors into the uncertainties in the rotation 
matrix by (22). The actual matrix of body direction vectors 
Bk+i and the actual attitude Ck+i also satisfy (23); 

Cl^iLk+i - Ll^iCk+i = 0, (27) 


where Lk+i = Ek+iWk+iBj.j^i G Using the iden¬ 

tity, S'(a;)Gl + ATs'(a;) = 5'({tr [A]/ 3 X 3 — A} x) for A G 
R^^^, X G R^, (27) can be written in the vector form 



i=l 


Then, we obtain 

m 

= (28) 

i=l 

where 

= - |tr [(cr+y L.+z] - (cr+y L.+zj 
w, {tr [y(e*)Tcr+,] / 3 X 3 - (29) 

This equation expresses the error in the measured attitude as 
a linear combination of the directional sensor errors. 


The perturbation of the angular velocity is equal 

to the angular velocity measurement error Vk+i- Substituting 
(26) into (13), we obtain 

5u;^+i=vk+i. (30) 

Define the error states xt+i = [(cr+/)^- e 

R®. Using (28) and (30), 

m 

xT+i=H,J2ATM^i+H,Vk+i, 

2=1 


where iJi — [I 3 X 3 , 03 x 3]"^;-^^2 — [ 03 x 3 ,-^ 3 x 3 ]"^ G R®^^ 
which expresses as a linear combination of the sensor 
errors A and v. From (15) and (16), each term on the right 
hand side is in the following uncertainty ellipsoids: 


HiAlXA+i e ^Q,H,AlXSUi [AXX) 

H2Vk+i G £r6 H2Tk+iHj^ . 


The measurement update finds a minimal ellipsoid containing 
the vector sum of these uncertainty ellipsoids. Expressions 
for a minimal ellipsoid containing multiple ellipsoids are 
given in [17] and [18], and PJ^i is given by 




In summary, the measured uncertainty ellipsoid at the k + 
(th step is defined by (23), (24), and (31); 

{Ck+uUk+i) G S{dX+i,^X+uPX+i)- (32) 

C. Filtering procedure 

The filtering procedure is to find a new uncertainty el¬ 
lipsoid compatible with the predicted and the measured 
uncertainty ellipsoids. From (21) and (32), the state at k + kh 
step lies in the intersection 

{Ck+l,UJk+l) G 

£{CX+^,u:X+i,Pr+i)- (33) 

Since it is inefficient to describe an irregular subset like the 
intersection of two ellipsoids in the state space numerically, 
we find a minimal uncertainty ellipsoid containing the inter¬ 
section. We omit the subscript {k + 1) in this subsection for 
convenience. 



The measurement uncertainty ellipsoid, ti)’”, P™), 

is identified by its center and the uncertainty 

ellipsoid in R®: 

(r,^a;™)e£:Re(06xi,P’”), (34) 

where = logm S so(3), ffw™ = lo - 

oj™ £ R3. Similarly, the predicted uncertainty ellipsoid, 
,u}^, P^), is identified by its center {C -^and the 
uncertainty ellipsoid in R®: 

iC^,SLo^)G£M<06xi,P^'), (35) 

where S'(C'^) = logm G so(3), Soj-I' = uj — u;f G 

R3. 

Define ("^^,6(1}'^^ G R^ such that 

Qf (36) 

a,/+ (37) 

Thus, gives the difference between the centers 

of the two ellipsoids. Using (36) and (37) we get 

67/=c-eS(c“UeS(CU^ 

(38) 

w/ = + {dto^f + div^) , (39) 

where we assumed that sufficiently small. 

Thus, the uncertainty ellipsoid obtained by the flow update, 
£{C^,pf) is given by the center (C'’",a)’”) of the 
measurement uncertainty ellipsoid and 

(40) 

r T I'T 

where S R®. 

We seek a minimal ellipsoid that contains the intersection 
of two uncertainty ellipsoids in R®: 

C ^Ra(i,P), (41) 

where x = (5(1)'^]'^ G R®. We obtain x and P as 

X = Lx'^f, 

P = P{q)iI-L)P^, 

where 

P{q) = l + q- (x’^ff{P"^)-^Lx^f, 

L = P^{P^ + q-^pf)-\ 

The constant q is chosen such that tr [P] is minimized. We 
convert x to points in TSO(3) using the common center 

In summary, a new uncertainty ellipsoid at the k + Zth step 
is defined by 

{Ck+i,0Jk+i) & £{Ck+i,tOk+i, Pk+i), (42) 


where 


Ck+i = cr+;e^(«), 

(43) 

LOk+l = 

(44) 

Pk+l = P- 

(45) 


The entire estimation procedure is repeated. The new un¬ 
certainty ellipsoid is used to predict the uncertainty ellipsoid 
till the next measurements are available, and the measure¬ 
ment update and the filtering procedures are performed. 
The center of the new uncertainty ellipsoid provides point 
estimates of the attitude and the angular velocity at the 
k + /th step. The uncertainty matrix represents the ellipsoidal 
bound on uncertainty. The size of the uncertainty matrix 
characterizes the accuracy of the estimates. If the size is 
small, we conclude that the estimates are accurate. This 
estimation scheme is optimal since the size of the new 
uncertainty ellipsoid is minimized. The eigenvector of the 
uncertainty matrix corresponding to the maximum eigenvalue 
shows the direction of the maximum uncertainty. 

V. Numerical Simulation 

Numerical simulation results are given for the estimation 
of the attitude dynamics of an uncontrolled rigid spacecraft in 
a circular orbit about a large central body, including gravity 
gradient effects. The detailed description of the on orbit 
spacecraft model is presented in [19]. 

The inertia of the spacecraft is chosen as J = 
diag [1, 2.8, 2], where overlines denote normalized variables. 
The maneuver is an arbitrary large attitude change completed 
in a quarter of the orbit, Tf = ^ s. The initial conditions 
are chosen as 

Co = diag[-l, -1, 1], wo = [2.3160, 0.4468, -0.5910]’^, 
Co = hx3, ^0 = [2.1160, 0.5468, -0.8910]’^. 

The corresponding initial estimation errors are ||(^o|| = 
180 deg, ||()wo|| = 21.43 rad/s. The initial uncertainty 
matrix is given by 

Po = 2diag (l80^)'[l, 1, 1], (30^)'[1, 1, 1] , 

so that XqPq^xq = 0.7553 < 1. 

We assume that measurements are available ten times in 
a quarter orbit. The measurement uncertainty matrices are 
given by 

SI = (^^io) -^3X3 rad^, Tk = (7^) /sxsradVs^ 

Fig. 2 shows simulation results for a typical realization of the 
bounded uncertainties, where the plot on the left shows the 
attitude estimation error and the angular velocity estimation 
error, and the right plot shows the size of the uncertainty 
ellipsoid. The estimation errors and the size of uncertainty 
decrease fast after the first estimation. The terminal attitude 
error is less than 1 deg. 
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Fig. 2. Attitude and angular velocity estimation errors with measurements 


VI. Conclusions 

The attitude estimation scheme presented here has no 
singularities since the attitude is represented by a rotation 
matrix, and the structure of the group of rotation matrices 
is preserved since it is updated by group operations in 
SO (3) using the Lie group variational integrator. The attitude 
estimator is also robust to the distribution of the uncertainty 
and the sensor noise, since it is based on deterministic 
ellipsoidal bounds on the uncertainty. The effects of process 
noise can be included by modifying the prediction procedure. 

Although not presented in this paper, we have obtained 
results for the modification of this scheme to the case when 
angular velocity measurements are not available. We intend 
to extend this estimation scheme to the combined attitude 
control and estimation problem for a rigid body in an attitude 
dependent potential, with the inclusion of process noise or 
disturbance forces. These topics will be dealt with in a future 
journal paper. 
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